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Abstract. Anisotropic flows (v2 and V4) of hadrons and light nuclear clusters are studied by 
a partonic transport model and nucleonic transport model, respectively, in ultra-relativistic and 
intermediate energy heavy ion collisions. Both number-of-constituent-quark scaling of hadrons, 
especially for isf meson which is composed of strange quarks, and number-of-nucleon scaling of light 
nuclear clusters are discussed and explored for the elliptic flow (V2). The ratios of V4/V2 of hadrons 
and nuclear clusters are, respectively, calculated and they show different constant values which are 
independent of transverse momentum. The above phenomena can be understood, respectively, by 
the coalescence mechanism in quark-level or nucleon-level. 

INTRODUCTION 

Anisotropic flow is of interesting subject in theoretical and experimental investigations 
on heavy ion collision dynamics [1, 2, 3, 4, 5, 6, 7, 8, 9, 10]. Many studies of the 1-th 
and 2-nd anisotropic flows, namely the directed flow and elliptic flow, respectively, re- 
vealed much interesting physics about the properties and origin of the collective flows. 
In particular, recent ultra-relativistic Au -1- Au collision experiments demonstrated the 
number of constituent-quark (NCQ) scaling from the transverse momentum dependence 
of the elliptic flow for different mesons and baryons at the Relativistic Heavy Ion Col- 
lider (RHIC) in Brookhaven National Laboratory (BNL) [11], it indicates that the par- 
tonic degree of the freedom plays a dominant role in formation of the dense matter in 
the early stage of collisions. Several theoretical models were proposed to interpret the 
NCQ-scaling of hadrons at RHIC [12, 13, 14, 15, 16]. Of which, a popular interpre- 
tation is assuming that mesons and baryons are formed by the coalescence or recom- 
bination of the constituent quarks. While at intermediate energy heavy ion collisions, 
the coalescence mechanism has been also used to explain the formation of light frag- 
ments and their spectra of kinetic energy or momentum of light particles some times 
ago [17, 18, 19, 20, 21]. A few studies investigated the mass dependence of the directed 
flow [22, 23]. However, systematic theoretical studies on the elliptic flow of different nu- 
clear fragments in intermediate energy domain in terms of the coalescence mechanism 
is rare. 

In this work, we will investigate the flow scaling of hadrons and of nuclear fragments 
with different microscopic transport theories, namely the partronic transport model 
in ultra-relativistic energies and nucleonic transport model at intermediate energies, 
respectively. A common naive coalescence mechanism in quark level or nucleonic 



level was applied to analyze the flow. Hadrons, namely baryons and mesons, can be 

formed by the quark-level coalescence. While, nuclear fragments can be formed by the 
nucleon-level coalescence. Dependences of the number-of-constituent-quark or number- 
of-nucleon (NN) of the elliptic flow are surveyed, respectively, and the ratio of V4 and 
is studied. 

Two examples are given to illustrate the NCQ or NN scaling of the elliptic flow in 
ultra-relativistic and intermediate energy HIC, respectively. A Multi-Phase Transport 
(AMPT) model was applied to simulate Au-i-Au at = 200 GeV/c for the former case 
and Isospin-Dependent Molecular Dynamics (IDQMD) was used to simulate ^^Kr + 
^^"^Sn at Ei,gam/A = 25 MeV/u for the later case. 



DEFINITION OF ANISOTROPIC FLOW AND MODEL 

INTRODUCTION 



Definition of anisotropic flow 

Anisotropic flow is defined as the different n-th harmonic coefficient v„ of an az- 
imuthal Fourier expansion of the particle invariant distribution [2] 
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where (|) is the azimuthal angle between the transverse momentum of the particle and 
the reaction plane. Note that in the coordinate system the z-axis along the beam axis, 
and the impact parameter axis is labelled as x-axis. The first harmonic coefficient vi 

represents the directed flow, vi = {cos^) = (^), where pt = \J~^l-\-^ is transverse 

momentum. V2 represents the elliptic flow which characterizes the eccentricity of the 
particle distribution in momentum space, 
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V2 = W2(^)) = (^^), (2) 
and V4 represents the 4-th momentum anisotropy, 
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Partonic/hadronic transport model: AMPT model 

AMPT (A Multi-Phase Transport) model [24] is a hybrid model which consists of 
four main processes: the initial condition, partonic interactions, the conversion from 
partonic matter into hadronic matter and hadronic interactions. The initial condition, 
which includes the spatial and momentum distributions of minijet partons and soft 



string excitation, are obtained from the HIJING model [25]. Excitation of strings melt 

strings into partons. Scatterings among partons are modelled by Zhang's parton cascade 
model (ZPC) [26], which at present includes only two-body scattering with cross section 
obtained from the pQCD with screening mass. In the default version of the AMPT model 
[27] partons are recombined with their parent strings when they stop interaction, and the 
resulting strings are converted to hadrons by using a Lund string fragmentation model 
[28]. In the string melting version of the AMPT model (we briefly call it as "the melting 
AMPT" model)[29], a simple quark coalescence model is used to combine partons into 
hadrons. Dynamics of the subsequent hadronic matter is then described by A Relativistic 
Transport (ART) model [30]. Details of the AMPT model can be found in a recent 
review [24]. It has been shown that in previous studies [29] the partonic effect can not 
be neglected and the melting AMPT model is much more appropriate than the default 
AMPT model when the energy density is much higher than the critical density for the 
QCD phase transition [24, 29, 16]. In the present work, the partonic interacting cross 
section in AMPT model with string melting is selected as lOmb. 



Nucleonic transport model: QIVID IVIodel 

The Quantum Molecular Dynamics (QMD) approach is an n-body theory to describe 
heavy ion reactions from intermediate energy to 2 A GeV. It includes several important 
parts: the initialization of the target and the projectile nucleons, the propagation of nu- 
cleons in the effective potential, the collisions between the nucleons, the Pauli blocking 
effect and the numerical tests. A general review about QMD model can be found in 
[31]. The IDQMD model is based on QMD model affiliating the isospin factors, which 
includes the mean field, two-body nucleon-nucleon (NN) collisions and Pauli blocking 
[32, 33, 34, 35, 36]. 

In the QMD model each nucleon is represented by a Gaussian wave packet with a 
width \fL (here L = 2.16 fm^) centered around the mean position riit) and the mean 
momentum Piit), 
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The nucleons interact via nuclear mean field and nucleon-nucleon collision. The 
nuclear mean field can be parameterized by 



f/(p,x,)=a(f )+|3(f )^+i(l-x,)K + Q,«^^^^-^x, + f/^"* (5) 

Po Po 2 Po 



with Po the normal nuclear matter density (here, 0. 16 fm~^ is used), p, p„ and Pp are the 
total, neutron and proton densities, respectively. is zth component of the isospin degree 
of freedom, which equals 1 or -1 for neutrons or protons, respectively. The coefficients 
a, P and y are parameters for nuclear equation of state (EOS). Cgy^ is the symmetry 
energy strength due to the difference of neutron and proton. In the present work, we take 
a = -124 MeV, P = 70.5 MeV and y = 2.0 which corresponds to the so-called hard EOS 



with an incompressibility of ^ = 380 MeV and Cgym = 32 MeV [31]. Vc is the Coulomb 
potential, f/^"^ is Yukawa (surface) potential. 

The time evolution of the colliding system is given by the generalized variational 
principal. Since the QMD can naturally describe the fluctuation and correlation, we 
can study the nuclear clusters in the model [31, 32, 33, 34, 35, 36]. In QMD model, 
nuclear clusters are usually recognized by a simple coalescence model: i.e. nucleons are 
considered to be part of a cluster if in the end at least one other nucleon is closer than 
fmin < 3.5 fm in coordinate space and pmin < 300 MeV/c in momentum space [31]. This 
mechanism has been extensively applied in transport theory for the cluster formation. 



NUIVIBER-OF-CONSTITUENT-QUARK SCALING OF THE 
ELLIPTIC FLOW AT RHIC ENERGIES 

In non-central Au+Au collisions at RHIC energies, the overlap region is anisotropic 
(nearly elliptic). Large pressure built up in the collision center results in pressure gradi- 
ent which is dependent on azimuthal angle, which generates anisotropy in momentum 
space, namely elliptic flow. Once the spatial anisotropy disappears due to the anisotropic 
expansion, development of elliptic flow also ceases. This kind of self-quenching process 
happens quickly, therefore elliptic flow is primarily sensitive to the early stage equation 
of state (EOS) [37]. 

Since many particles have been investigated already for demonstrating NCQ scaling 
of elliptic flow at RHIC. However, the study of strangeness particle, ^, is rare. Figure 1 
shows both the measurement and the calculation of (|)-meson elliptic flow V2 at RHIC 
[38, 39, 40, 41]. For comparison, V2 of K^+K~ and p + p from [42] are plotted together. 
The most striking thing in the figure is that the (t)-meson has a non-zero V2 in the 
hydrodynamic region (px < 2 GeV/c) and it has a significant V2 signal for pj > 2 
GeV/c which is comparable to that of the K^+K . Since the formation of (|)-meson 
through kaon coalescence at RHIC has been ruled out by previous STAR measurement 
[43], and the low interaction cross-section of (t)-meson makes the contribution to flow 
due to hadronic re-scattering processes very small [24], it may be possible to directly 
measure the flow of strange quark via the flow of (|)-meson. Additionally, since the (|)- 
meson V2 values as a function of pr are similar to those of other particles, this means 
that the heavier 5-quarks flow as much as the lighter u and d quarks. This can happen if 
there are a significant number of interactions between the quarks before hadronization. 
It is the signal of partonic collectivity of the system! A prediction [16] from the naive 
coalescence model which is implemented in A Multi-Phase Transport (AMPT) model 
with the string melting scenario [44] is also plotted here. The model describes the data 
quite well at pj > 1-5 GeV/c. This is another hint of coalescence of strange quark as 
dominant production mechanism for (])-meson. Further evidence can be obtained from 
the Number-of-Constituent-Quark (NCQ) picture in Figure 1. 

The RHIC experimental data demonstrated a scaling relationship between 2-nd flow 

n/2 

(v2) and n-th flow (v„), namely Vn{pt) ~ (pt) [45]. It has been shown [46, 47] 
that such scaling relation follows from a naive quark coalescence model [13] that only 
allows quarks with equal momentum to form a hadron. Denoting the meson anisotropic 
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FIGURE 1. Left panel: pr dependence of V2 for (|)-meson which is compared with the cases of K^+K^ 
and p+'p- AMPT calculation are taken from Ref. [16]. Right panel: NCQ scaled V2 as a function of NCQ 
scaled transverse momentum. Lines represent NCQ-scaling parameterization. 



flows by Vn,M{pt) and baryon anisotropic flows by Vn^siPt), Kolb et al. found that 

V4,M(Pf) = (l/4)v2M(/'f) for mesons and V4,B{pt) = (1/3)v2.5(a) forbaryons if quarks 
have no higher-oraer anisotropic flows. Since mesons dominate the yield of charged 
particles in RHIC experimental data, the smaller scaling factor of 1 /4 than the empirical 
value of about 1 indicates that higher-order quark anisotropic flows cannot be neglected. 
Including the latter contribution, one can show that 

V4,M ^4,g I (f,\ 

v2 ~ 4 2 v2 ' v2 ~ 3 3 v2 ' 

where Vn,q denotes the quark anisotropic flows. The meson anisotropic flows thus satisfy 
the scaling relations if the quark anisotropic flows also satisfy such relations. Left panel 
of Fig. 2 shows the AMPT results for V2, V4, and vs from minimum bias events for 
Au+Au at ^/s = 200 A Gev with parton scattering cross section 10 mb [47]. It is 
seen that the parton anisotropic flows from the AMPT model indeed satisfy the scaling 

n /2 

relation Vn{pt) ~ (pt) [47]. However, as shown in the middle panel of Fig. 2, V4 is 
experimentally larger than the simulation and the ratio of V4/V2 is about 1.2 for charged 
hadrons [45, 48] (see right panel of Fig. 2), which means that the fourth-harmonic flow 
of quarks V4 must be greater than zero. One can go one step further and assume that the 
observed scaling of the hadronic V2 actually results from a similar scaling occurring at 
the partonic level. In this case, if one assumes [46, 47] that the scaling relation for the 
partons is as follows: 

vi = (v^)^ (7) 

and then hadronic ratio V4 / then equals 1/4+1/2 = 3/4 for mesons and 1/3 + 1/3 = 
2/3 for baryons, respectively. However, since this value is measured to be 1.2 for 
charged hadrons and pions, so that either the parton scaling relation [Eq. (7)] must 
have a proportionality constant of about 2, or the simple coalescence model needs 
improvement. 




FIGURE 2. (Color online) Left panel (AMPT calculations): Transverse momentum dependence of 
midrapidity parton anisotropic flows V2, V4, and vg from minimum bias events for Au+Au at = 200 A 
Gev with parton scattering cross section 10 mb. Also plotted are (open triangles) and Vj (open stars). 
Figure taken from [47]; Middle panel (STAR data): Minimum bias measurements of anisotropic flow of 
charged hadrons for different harmonics for Au+Au at ^/s = 200 A GeV. The dashed lines show 1 .2 
and 1.2v2, respectively. Figure taken from [45]; Right panel (STAR data): Ratio of V4/V2 vs pr of charged 
hadrons for Au+Au at ^ = 200 A GeV . The straight line represents the mean of all entries at a value of 
1.2 [48]. 



NUMBER-OF-NUCLEON SCALING OF THE ELLIPTIC FLOW 
AT INTERMEDIATE ENERGIES 

In intermediate energy, we simulated ^^Kr + ^^^Sn at 25 IVIeV/nucleon and impact pa- 
rameter of 7 - 10 fm with IQIVID [49] . 50,000 events have been accumulated. The system 
tends to freeze-out around 120 fm/c. In this work, we extract the following physics re- 
sults at 200 fm/c. Fig. 3(a) shows transverse momentum dependence of elliptic flows for 
mid-rapidity light fragments. From the figure, it shows that the elliptic flow is positive 
and it increases with the increasing pt. It reflects that the light clusters are preferentially 
emitted within the reaction plane, and particles with higher transverse momentum tend 
to be strongly emitted within in-plane , i.e. stronger positive elliptic flow. In compari- 
son to the elliptic flow at RHIC energies, the apparent behavior of elliptic flow versus 
Pt looks similar, but the mechanism is obviously different. In intermediate energy do- 
main, collective rotation is one of the main mechanisms to induce the positive elliptic 
flow [50, 51, 52, 53, 54, 55]. In this case, the elliptic flow is mainly driven by the at- 
tractive mean field. However, the strong pressure which is built in early initial geomet- 
rical almond-type anisotropy due to the overlap zone between both colliding nuclei in 
coordinate space will rapidly transforms into the azimuthal anisotropy in momentum 
space at RHIC energies [11]. In other words, the elliptic flow is mainly driven by the 
stronger outward pressure. Fig. 3(b) displays the elliptic flow per nucleon as a func- 
tion of transverse momentum per nucleon, and it looks that there exists the number of 
nucleon scaling when pt/A < 0.25 GeV/c. This behavior is apparently similar to the 
number of constituent quarks scaling of elliptic flow versus transverse momentum per 
constituent quark (pt/n) for mesons and baryons which was observed at RHIC [11] and 
shown in Fig. 3(b). 

Recognizing the V4/V2 behaviors at RHIC energies, we would like to know what 
the higher order momentum anisotropy in the intermediate energy is. So far, there is 
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FIGURE 3. (a) Elliptic flow as a function of transverse momentum {pt). Squares represent for neutrons, 
circles for protons, triangles for fragments of A = 2, diamonds for A = 3 and stars for A = 4; (b) Elliptic 
flow per nucleon as a function of transverse momentum per nucleon. The symbols are the same as (a). 

neither experimental data nor theoretical investigation for the higher order flow, such 
as V4, in this energy domain. In the present work, we explore the behavior of V4 in 
the model calculation. Fig. 4 shows the feature of V4. Similar to the relationship of 
V2/A versus pt/A, we plot V4/A as a function of pt/A. The divergence of the different 
curves between different particles in Fig. 4(a) indicates no simple scaling of nucleon 
number for 4-th momentum anisotropy. However, if we plot V4/A^ versus (pt/A)^, it 
looks that the points of different particles nearly merge together and it means a certain 
of scaling holds between two variables. Due to a nearly constant value of V4/V2 in the 
studied pt range (see Fig. 4(c)) together with the number-of-nucleon scaling behavior of 
V2/A vs Pt/A, V4/A^ should scale with (pt/A)^, as shown in Fig. 4(b). If we assume 
the scaling laws of mesons (Eq. 6) are also valid for A = 2 and 3 nuclear clusters, 
respectively, then V4/V2 for A = 2 and 3 clusters indeed give the same value of 1/2 as 
nucleons, as shown in Fig. 4(c). Coincidentally the predicted value of the ratio of V4/V2 
for hadrons is also 1/2 if the matter produced in ultra-relativistic heavy ion collisions 
reaches to thermal equilibrium and its subsequent evolution follows the laws of ideal 
fluid dynamics [56]. It is interesting to note the same ratio was predicted in two different 
models at very different energies, which is of course worth to be further investigated 
in near future. Overall speaking, we learn that V4/V2 is approximately 1/2 in nucleonic 
level coalescence mechanism, which is different from 3/4 for mesons or 2/3 for baryons 
in partonic level coalescence mechanism. 

SUMMARY 

In summary, we applied AMPT and IDQMD model to investigate the behavior of 
anisotropic flows, namely V2 and V4, versus transverse momentum for hadrons in Au+Au 
at y/s^ = 200 GeV/c as well as light nuclear fragments for 25 MeV/nucleon ^^Kr + ^^^Sn 
collisions at 7- 10 fm of impact parameters. Both V2 and V4 generally show positive values 
and increase with p; in two very different energies. By the number scaling of constituent- 
quarks or nucleons, the curves of elliptic flow for different hadrons or nuclear fragments 
approximately collapse on the similar curve, respectively, which means that there exists 
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FIGURE 4. (a) V4/A as a function of 75, /A for different particles, namely, neutrons (squares), protons 
(circles), fragments of A = 2 (triangles), A = 3 (diamonds) and A = 4 (stars), (b) v^jA? as a function of 
{pt /A)^. (c) the ratios of vi\,jv\ for different particles vs pt- 



an elliptic flow scaling on the constituent quark number or nucleon number. NCQ 
scaling stems from the partonic coalescence and NN scaling originates from nucleonic 
coalescence. For 4-th momentum anisotropy V4, it seems to be scaled by v\, and V4/V2 
is different for hadrons and nuclear fragments. AMPT predicts a smaller value of va,Iv\ 
than the data ~ 1.2 and the value of nuclear fragments ~ 0.5. The former indicates that 
either the parton scaling relation [Eq. (7)] must have a proportionality constant of about 
2, or the simple quark coalescence model needs improvement. For the later, the data is 
unable to be available yet. Therefore it will be of very interesting if one can measure this 
ratio in intermediate energy HIC. 
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